d/u Asymmetry and the Origin of the Nucleon Sea 



OO 
On 

U . 

< 



> 
OO 
00 
(N 

^- 

o 

00 
Oh! 



X 



J.C. Peng^, G.T. Garvey/, T.C. Awes 1 , M.E. Beddo' 1 , M.L. Brooks/, C.N. Brown c , J.D. Bush a , T.A. Carey/, 
T.H. Chang' 1 , W.E. Cooper c , C.A. GagliardP, D.F. Geesaman 6 , EA. Hawker 3 ', X.C. He d , L.D. Isenhower , 
S.B. Kaufman 6 , D.M. Kaplan 6 , P.N. Kirk^, D.D. Koetke fc , G. Kyle 71 , D.M. Lee/, W.M. Lee d , M.J. Leitch/, 
N. Makinsfl, P.L. McGaughey/ , J.M. Moss/, B.A. Mueller 6 , P.M. Nord fc , B.K. Park/, V. Papavassiliou' 1 , G. Petitt d , 
P.E. Reimer/, M.E. Sadler J. Selden' 1 , P.W. StanW, W.E. Sondheim/, T.N. Thompson/, R.S. TowelpQ, 

R.E. Tribble 3 , MA. Vasiliev 3 !, Y.C. Wangf, Z.F. Wang°, J.C. Webb' 1 , J.L. Willis , D.K. Wise a , G.R. Young* 

(FNAL E866/NuSea Collaboration) 

a Abilene Christian University, Abilene, TX 79699 
b Argonne National Laboratory, Argonne, IL 60439 
c Fermi National Accelerator Laboratory, Batavia, IL 60510 
d Georgia State University, Atlanta, GA 30303 
e Illinois Institute of Technology, Chicago, IL 60616 
f Los Alamos National Laboratory, Los Alamos, NM 87545 
9 Louisiana State University, Baton Rouge, LA 70803 
h New Mexico State University, Las Cruces, NM, 88003 
1 Oak Ridge National Laboratory, Oak Ridge, TN 37831 
3 Texas A & M University, College Station, TX 77843 
k Valparaiso University, Valparaiso, IN 46383 
(February 1, 2008) 

The Drell-Yan cross section ratios, a(p + d)/a(p + p), measured in Fermilab E866, have led to 
the first determination of d(x)/u(x), d(x) — u(x), and the integral of d(x) — u(x) for the proton 
over the range 0.02 < x < 0.345. The E866 results are compared with predictions based on parton 
distribution functions and various theoretical models. The relationship between the E866 results 
and the NMC measurement of the Gottfried integral is discussed. The agreement between the E866 
results and models employing virtual mesons indicates these non-perturbative processes play an 
important role in the origin of the d, u asymmetry in the nucleon sea. 
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Recent measurements Jl)-^ have revealed a marked 
asymmetry in the distributions of up and down quarks in 
the nucleon sea. While no known symmetry requires u to 
equal d, a large d/u asymmetry was not anticipated. The 
principal reason for expecting symmetry between up and 
down quarks in the sea is an assumption that the sea orig- 
inates primarily from q-q pairs produced from gluons. As 
the masses of the up and down quarks are small compared 
to the confinement scale, nearly equal numbers of up and 
down pairs should result. Indeed, a theoretical investi- 
gation H| of the light-quark asymmetry in the nucleon 
concluded that perturbative processes do not give rise to 
asymmetries in the up, down sea exceeding 1%. Thus 
a large d/u asymmetry requires a non-perturbative ori- 
gin for an appreciable fraction of these light antiquarks. 
This paper draws together several implications arising 
from this observed d/u asymmetry - the effect of these 
measurements on existing parton distributions, an exam- 
ination of the compatibility of the measurements of this 
asymmetry, and the origin of the effect. 

The issue of the equality of u and d was first encoun- 
tered in measurements of the Gottfried integral ||, de- 
fined as 



[FZ(x,Q 2 )-F?(x,Q 2 )] /xdx, 



(1) 



where Ff and F£ are the proton and neutron structure 
functions measured in deep inelastic scattering (DIS) ex- 
periments. Iq can be expressed in terms of the valence 
and sea quark distributions of the proton as: 



Ic 



[u v (x, Q 2 ) - d v (x, Q 2 )] dx 

+ | [u(x, Q 2 ) - d(x, Q 2 )} dx. (2) 

Under the assumption of a u, d flavor-symmetric sea in 
the nucleon, the Gottfried Sum Rule (GSR) @, I G = 1/3, 
is obtained. Measurements of muon DIS on hydrogen and 
deuterium by the New Muon Collaboration (NMC) deter- 
mined that f ° 8 04 [F£(x) ~ F%{x)\ /x dx = 0.221 ± 0.021 
at Q 2 = 4 GeV 2 M. Extrapolating to x = through 
the unmeasured small-x region, the Gottfried integral is 
projected to be 0.235 ± 0.026, significantly below 1/3. 

Although the violation of the GSR observed by NMC 
can be explained by assuming pathological behavior of 
the parton distributions at x < 0.004, a more natural 
explanation is to abandon the assumption u = d. Specif- 
ically, the NMC result implies 
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FIG. 1. Comparison of the E866 d — u results with the 
predictions of various models as described in the text. 



[d{x) - u(x)] dx = 0.148 ± 0.039. 



(3) 



The Fermilab E866 measurement |J of the ratio of Drell- 
Yan yields from hydrogen and deuterium directly de- 
termines the ratio d(x)/u{x) for 0.02 < x < 0.345. An 
excess of d over u is found over this x range, supporting 
the observation by NMC that the GSR is violated. 

The d/u ratios measured in E866, together with the 
CTEQ4M § values for d + u, were used to obtain d — u 
over the region 0.02 < x < 0.345 (Fig. |l|). As a flavor 
non-singlet quantity, d{x) — u(x) has the property that 
its integral is Q 2 -independent |8J. Furthermore, it is a di- 
rect measure of the contribution from non-perturbative 
processes, since perturbative processes cannot cause a 
significant d, u difference. As shown in Fig. 1, the x de- 
pendence of d — u at Q — 7.35 GeV can be approximately 
parametrized as 0.05ar a5 (l - x) 14 (l + lOOx). 

Integrating d(x) — u(x) from E866, one finds 
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[d(x) — u(x)] dx = 

0.068 ± 0.007(stat.) ± 0.008(syst.) (4) 



at Q — 7.35 GeV ||. To investigate the compatibility 
of this result with the NMC measurement (Eq. 3), the 
contributions to the integral from the regions x < 0.02 
and x > 0.345 must be estimated. Table | lists the val- 
ues for the integral of d — u over the three regions of 
x for three different parton distribution function (PDF) 
parametrizations at Q = 7.35 GeV. For x > 0.345, the 
contribution to the integral is small (less than 2%). The 
three parametrizations predict that the bulk of the con- 
tribution to the integral comes from 0.02 < x < 0.345. 
Since CTEQ4M provides a reasonable description of the 
E866 data in the low- a; region ||], and the contribution 
from the high-cc region is small, we have used CTEQ4M 
to estimate the contributions to the integral from the un- 
measured x regions. This procedure results in a value 



TABLE I. Values for J [d(x] - u(x)] dx over various x 
ranges, evaluated at Q = 7.35 GeV, for various PDF 
parametrizations. Values deduced from E866 are also listed. 



x range CTEQ4M MRS(R2) GRV94 



E866 



0.345 - 1.0 
0.02 - 0.345 
0.0 - 0.02 
0.0 - 1.0 



0.00192 
0.0765 
0.0296 
0.1080 



0.00137 
0.1011 
0.0588 
0.1612 



0.00148 
0.1027 
0.0584 
0.1625 



0.068 ± 0.011 



0.100 ± 0.018 



/o [d(x) - u(x)] dx = 0.100 ± 0.007 ± 0.017, which is 
2/3 the value deduced by NMC. The systematic error in- 
cludes the uncertainty (±0.015) due to the unmeasured x 
regions, estimated from the variation between CTEQ4M 
and MRS(R2) (jlj. This result is consistent with the 
integral of the parametrized fit shown in Fig. 1. 

The difference between the NMC and E866 results for 
the d — u integral raises the question of the compatibil- 
ity of the two measurements. Figure || shows the NMC 
data for F% — F% at Q = 2 GeV, together with the fits 
of MRS(R2) and CTEQ4M. Both PDF parametrizations 
give very similar results for F% — F% . However, their 
agreement with the NMC data is poor, especially in the 
region 0.15 < x < 0.4. It is instructive to decompose 
F£( x ) ~ F-2( x ) m t° contributions from valence and sea 
quarks: 



\x \u v (x) 



d v {x)} 



fx [u{x) 



d(x)]. (5) 



Two PDF parametrizations of these contributions are 
shown in Fig. ^|. The valence contribution is positive, 
while the contribution from the sea is negative. The 
MRS(R2) and CTEQ4M parametrizations give notice- 
ably different values for the valence and sea contribu- 
tions, though their net results for F^ — F£ are very sim- 
ilar. As shown in Fig. [| the E866 data provide a direct 
determination of the sea-quark contribution to F% — K^, 
and can be used to distinguish between different PDF 
parametrizations that produce similar fits to the NMC 
data. As the direct determination of d{x)—u{x) is smaller 
than obtained from either PDF set, the parameters must 
be adjusted to reduce the magnitude of both the sea and 
valence distributions in the interval 0.03 < x < 0.3. This 
reduction will force an increase in the valence contribu- 
tion to the integral from x < 0.03 and could therefore 
bring the results from E866 and NMC into better accord. 
Fig. U also suggests that the reason for the difference be- 
tween the PDF fits and the NMC results in the interval 
0.15 < x < 0.4 is that the PDFs cannot accommodate 
the rapid variation in the asymmetry of the nucleon sea 
as a function of x revealed by E866. 

The E866 data also allow the first determination of the 
momentum fraction carried by the difference of d and u. 
We obtain J^ 45 x [d{x) - u(x)] dx = 0.0065 ± 0.0010 
at Q = 7.35 GeV. If CTEQ4M is used to estimate the 
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FIG. 2. Ff — F^ 1 as measured by NMC compared with pre- 
dictions based on the CTEQ4M and MRS(R2) parametriza- 
tions. Also shown are the E866 results for the sea-quark con- 
tribution to F% — F£ ■ For each prediction, the top (bottom) 
curve is the valence (sea) contribution and the middle curve 
is the sum of the two. 



contributions from the unmeasured x regions, one finds 
that /o x [d(x) - u(x)] dx = 0.0075±0.0011, roughly 3/4 
of the value obtained from the PDF parametrizations. 
Unlike the integral of d(x)— u(x), the momentum integral 
is Q 2 -dependent and decreases as Q 2 increases. 

We now turn to the origin of the d/u asymmetry. As 
early as 1982, Thomas [[ll] pointed out that the virtual 
pions that dress the proton will lead to an enhancement 
of d relative to u via the (non-perturbative) "Sullivan 
process." Sullivan |Q previously showed that in DIS 
these virtual mesons scale in the Bjorken limit and con- 
tribute to the nucleon structure function. Following the 
publication of the NMC result, many papers [13-20] have 
treated virtual mesons as the origin of the asymmetry in 
the up, down sea of the nucleon. 

Using the notion that the physical proton (p) may 
be expanded in a sum of products of its virtual meson- 
baryon (MB) states, one writes p = (1 — a)po + aMB, 
where a is the probability of the proton being in vir- 
tual states MB and po is a proton configuration with a 

that 



symmetric sea. It is easy to show 14 



[d(x, Q 2 ) - u(x, Q 2 )} dx = (2a - 6)/3 



(6) 



where a is the probability of the virtual state ttN and 
b the probability for ttA. These two configurations are 
the dominant intermediate MB states contributing to 
the asymmetry |l8|,|l^]. Further, most recent calcula- 
tions of the relative probability of these two configura- 
tions find a w 26 |ttq,^9[. Using the value for the in- 
tegral extracted from E866 and assuming a = 26 yields 
a = 2b = 0.20 ± 0.036, requiring a substantial presence 



of virtual mesons in the nucleon in this model. Follow- 
ing the observation [|l7j that these configurations have a 
large impact on the spin structure of the nucleon because 
pion emission induces spin flip, one can show that 

Au p - Ad p = 5/3 - 20(2a + b)/27 = 1.296 ± 0.067 (7) 

using the above values of a and b determined from the 
E866 result. Here Au p (Ad p ) is the total spin carried by 
up (down) quarks in the proton. This value is in good 
agreement with the measured axial coupling constant for 
the nucleon, gA — 1.260 ± 0.003 plj], and increases con- 
fidence in the virtual meson-baryon picture. 

The x dependences of d — u and d/u obtained in 
E866 provide important constraints for theoretical mod- 
els. Fig. [I] compares d{x) — u(x) from E866 with a virtual- 
pion model calculation following the procedure detailed 
by Kumano |Q . Since the E866 results are shown at a Q 
of 7.35 GeV, the SMRS(P2) |§2j parametrization for the 
pion structure functions at this Q is employed. The curve 
labeled "virtual pion A" in Fig. uses a dipole form with 
A = 1.0 GeV for the itNN and ttA^A form factors, and 
is seen to underpredict the magnitude of d — u. However 
as has been noted Jl8|,[l9| , A production experiments p3] 
suggest a considerably softer form factor for 7rA^ A than 
for ttNN. Indeed much better agreement with the E866 
results is obtained by reducing A for the irN A form fac- 
tor to 0.8 GeV, as shown by the curve labeled "virtual 
pion B" in Fig. jl]. This fit produces a value of 0.11 for 
the integral of d — u and 1.18 for gA- If A is chosen to 
be 0.9 GeV (0.7 GeV) for the ttAW \ttNA) form factor, 
one finds nearly exact accord with the values cited in the 
previous paragraph. 

A different approach for including the effects of vir- 
tual mesons has been presented by Eichten, Hinchliffe, 
and Quigg ||r^| and further investigated by Szczurek et 
al. poll . In the framework of chiral perturbation theory, 
the relevant degrees of freedom are constituent quarks, 
gluons, and Goldstone bosons. In this model, a portion of 
the sea comes from the couplings of Goldstone bosons to 
the constituent quarks, such as u — > dn + and d — * utt~ . 
The excess of d over u is then simply due to the additional 
up valence quark in the proton. The predicted d — u from 
the chiral model is shown in Fig. |] as the dotted curve. 
We follow the formulation of Szczurek et al. [^0) to cal- 
culate d(x) — u(x) at Q = 0.5 GeV, and then evolve the 
results to Q = 7.35 GeV. In the chiral model, the mean- 
x of d — u is considerably lower than in the virtual-pion 
model just considered. This difference reflects the fact 
that the pions are softer in the chiral model, since they 
are coupled to constituent quarks which on average carry 
only 1/3 of the nucleon momentum. The x dependence 
of the E866 data favors the virtual-pion model over the 
chiral model, suggesting that correlations between the 
chiral constituents should be taken into account. 

Another non-perturbative process that can produce a 
d, u asymmetry is the coupling of instantons to the va- 
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agreement between the E866 d — u data and the virtual- 
pion model indicates that virtual meson-baryon compo- 
nents play an important role in determining non-singlet 
structure functions of the nucleon. Future experiments 
extending the measurements of d/u to other x and Q 2 
regions can further illuminate the interplay between the 
perturbative and non-perturbative elements of the nu- 
cleon sea. 

This work was supported in part by the U.S. Depart- 
ment of Energy. 
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FIG. 3. Comparison of E866 d/u ratios at Q = 7.35 GeV 
with the predictions of two virtual pion models and the chiral 
model. The fit 1 + 1120a; 2 ' 75 (1 - a;) 15 is also shown. 



lence quarks. An earlier publication j24j] presented an 
asymmetry due to instantons but parametrized the re- 
sult in terms of the asymmetry observed in NMC, and 
therefore has no independent predictive power. Also the 
ad hoc x dependence used for d{x)/u(x) is in poor agree- 
ment with the E866 result. We are unable to determine 
if better agreement with data can be obtained by an im- 
proved parametrization within the instanton model. 

It is also instructive to compare the model predic- 
tions of d(x)/u(x) with the E866 results. Figure || shows 
that the two virtual-pion models and the chiral model 
give d(x)/u(x) values very different from the E866 re- 
sult. Note that these calculations do not include the 
perturbative processes g — > uu, dd which generate a sym- 
metric sea. Indeed, the d/u data provide valuable in- 
formation on the relative importance of the perturbative 
(symmetric) versus the non-perturbative sea. The chiral 
model predicts a d/u ratio of 11/7 for all x. Figure || 
shows that the observed ratio nearly equals this value 
for 0.1 < x < 0.2, leaving little room for any pertur- 
batively generated symmetric sea in this interval, which 
seems unreasonable. The same problem also arises for the 
virtual-pion model A. In contrast, the virtual-pion model 
B readily accommodates contributions from a symmetric 
sea. 

In summary, E866 has provided the first determination 
of d/u, d — u, and the integral of d — u over the range 
0.02 < x < 0.345. It provides an independent confirma- 
tion of the violation of the Gottfried Sum Rule reported 
from DIS experiments. The magnitude of the integral of 
d — u over the region 0.02 < x < 0.345 is smaller than 
obtained from some current PDF parametrizations. This 
indicates that the violation of the Gottfried Sum Rule is 
likely smaller than reported by NMC. Together with the 
NMC data, the E866 results impose stringent constraints 
on both sea- and valence-quark distributions. The good 
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